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Abstract 

The fluctuation dynamics in steady state fragmentation and aggregation processes is analyzed 
from a kinetic point of view. The fluctuations An(x,t) in the number of particles of size between 
x and x + Ax, in a sample with average occupation number AN(x), appear to behave like sums of 
uncorrelated events: < An(x,t) 2 >oc AJV(a;). The implications of this fact for the applicability of 
a mean field description, and possible relation with intermittency phenomena in three dimensional 
turbulence are discussed. 
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The analysis of fragmentation and agglomeration processes focuses usually on the determina- 
tion of the distribution N(x,t) for the particle size x, which gives the mean field approximation 
for the problem. Once an explicit expression for the transition probability a for the process is 
available, the distribution N can be obtained from a Boltzmann-like equation, which is linear for 
binary fragmentation (each break-up event results in only two fragments) [1-2]: 

/>GO i>X 

d t N(x,t) = 2 dya(y -+ x,y - x)N(y,t) - N(x,t) dya(x -+ y, x - y), (1) 

Jx Jo 

while, for agglomeration from binary collisions, one has the Smoluchowski equation [3]: 

1 r - 
d t N(x,t) = - J dya(y, x - y -> x)N(y, t)N(x - y, t) 

-N(x,t) dya(x,y ^ x + y)N(y,t). (2) 
Jo 

Clearly, while Eqn. (1) is exact, Eqn. (2) is an approximation, which neglects correlations and is 
based usually on assumptions of diluteness [4]. 

The field of application of Eqns. (1-2), is clearly very wide (see for instance [5] and [6] and 
references therein). The possibility of using a mean field description like the one provided by Eqns. 
(1) and (2) depends on whether fluctuations can be considered small. In the case of agglomeration, 
this affects the validity of the approximations leading to Eqn. (2), since large fluctuations mean 
that correlations in the form: < N(x,t)N(y,t) >, with < > indicating ensemble average, can 
dominate over the factorized terms N(x,t)N(y,t) =< N(x,t) >< N(y,t) >. 

The role of fluctuations in fragmentation was examined recently By Esipov et AL, in the case 
of break-up probabilities in the form: a(x y, x — y) = x^ [7]. Their analysis showed a situation 
of dominance of the fluctuations down the cascade in the shattering range [6], corresponding to 
fragmentation accelerating down the cascade in such a way that an initial object disintegrates into 
zero size fragments in a finite time. The relevance to intermittency issues in turbulence modeling 
(see e.g. [8]) is obvious. However, the reason why a transition to dominant fluctuations occurs in 
this way is not yet clear. 

The purpose of this letter is to clarify this issue, at least in situations of constant mass current, 
and to extend the analysis of [7] to agglomeration processes. Despite the differences between 
fragmentation and agglomeration processes, it will appear that the mechanism for the transition to 
a situation of fluctuation dominance, is the same in both. While in [7], the approach was essentially 
a statistical mechanical one, based on an analysis in "T-space" of the probability functional PfiV], 
the approach here will be to study directly the "//-space" kinetic equations for the N(x,t), whose 
mean field approximation is given by Eqns. (1) and (2). 

Fluctuations are analyzed identifying a set of mother particles, and defining a sample, as the 
subset of fragments (or agglomerates), which originate from them. The mother particles can be 
characterized for instance as the largest (in the case of fragmentation) or smallest (in the case of 
agglomeration) particles, which are present at any time, in a prescribed volume element V of the 
system. The volume V is chosen large enough to guarantee a sufficient sample, which, in the case 
of coalescence, is necessary also to minimize interactions with particles coming from the outside of 
V. In this way it will be possible also to consider n(x,t) = N(x,t) — N(x) as a spatial fluctuation. 

Binary fragmentation 

The kinetic equation for a fragmentation process can be written in the form, after discretizing 
in time and particle size: 

OO X — 1 

N(x,t+1) = N(x,t) + 2 W(y,x;t)-2^2w iX ,y;t), (3) 

y=x+l y=l 
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where W(x,y;t) = Sj=i' w ( x —> y, x — y,t;j), with w(x —> y, x — y,t;j) a stochastic variable 
that is 1 or 0, depending on whether the j-th particle of size x undergoes or not the breaking 
process: x y, x — y. If the break-up events occur independently, the probability that the sums 
W(x,y;t) take certain values, will be that of distributing W(x,y;t) particles, out of N(x,t), in 
slots y = 1, ...x — 1, with: W(x,y;t) = W(x,x — y;t). This is a multinomial distribution that, for 
N(x) 1, is approximated by: 

x y — 1 

where AW(x,y;t) = W(x,y;t) — a(x,y)N(x), x' = -| for x even and x' = for odd x, while 
a(x,y) = a(x — ► y, x — y) is the probability of an individual breaking event. Averaging Eqn. (3) 
over the distribution given by Eqn. (4) and taking the continuum limit, gives immediately Eqn. 
(1), while the fluctuations obey the stochastic (linear) integral equation: 

/>GO t>X 

d t n(x,t)-2 dya(y,x)n(y,t) + n(x,t) dya(x,y) = £(x,t) (5) 
Jx Jo 

where: 

/>oo 1>X 

e(M) = 2/ dyAW(y,x;t)- dyAW(x,y;t), (6) 

Jx Jo 

plays the role of external noise term. Focusing on the case of homogeneous breaking probabilities: 
a(x, y) = x@ f(y/x) with lim^o z f( z ) = 0, one obtains a stationary solution of Eqn. (1), requiring 
the mass current J(x) = 2 J dy J Q dza(y, z)zN(y) to be ^-independent. From power counting: 
N(x) ~ a; -7 , with 7 = 3 + (3. 

Given the statistics for AW, the correlation for the external noise £ can be calculated and one 
finds that it contains a singular piece: < £(s, t)£(x', t') > s — S(t — t')6(x — x')x~ 2 . The correlation 
for n, therefore, receives itself a singular contribution: < n(x,t)n(x' , /') > s , whose dynamics is 
dominated by the balance between < £(s, t)£(x', t') > s , and the loss term n(x,t) J dya(x,y) in 
eqn. (5). One obtains therefore: 

< n(x,t)n(x',t') > s cx 8(x - x')N(x) exp(-T(x)\t - t'\), (7) 

with r(a;) = x^ +1 dzf(z). This result implies that the fluctuations An(x,t) in the number of 
particles, in a small interval of size Ax, behave like the fluctuations in a sample of statistically 
independent events: < An(x) 2 >~ AN(x). From Eqn. (7), fluctuations become dominant down 

the cascade: < ^^^2 > 00 as a; ^ 0, when AN(x —> 0) = 0, i.e. for fi < —2; this coincides with 
the result obtained by Esipov et Al. [7] in the case of a Yule-Ferrin cascade [9]. 

These results are contingent upon neglecting any correlations among the break-up events. In 
particular, if one admits a finite correlation time in the fragmentation rate fluctuation t^(x,t), (which 
is especially important for turbulence modeling) different mechanisms for fluctuations generation 
become possible. The issue is already apparent, in as simple a turbulence model as a Yule-Ferrin 
cascade, in which each eddy splits into two equal fragments, conserving eddy energy and volume. 
From Kolmogorov scaling [10], one has for the energy of the eddies, at the ra-th step in the cascade: 

E m = N m x m ~ /m = 2 s m , where N m = AN{x m ), and x m = 2 m is the energy of the individual 
eddy in the ra-th generation, with linear size l m = 2~ ! 3 . This fixes also the break-up probability: 
a m = a(x m ,x m+ i) = E^ 1 . If, following [11], the finite correlation time of the break-up process is 
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kept into account, the amplitude of the external noise in Eqn. (5) is corrected by a factor the 
result: < An(x) 2 >~ N(x) ceases to be valid, and one obtains the equation for (f> m = iV~ 1 ra m : 



d t (j>m + a m d m c/) m = a m d m (a m 2 ^ 



where d m is a finite difference and < £ m (i)£ m '(i') >= S mm i S(t — t'). Solution of this equation gives 

then the result: < (f) 2 m >~ 1 + m(d m log a m ) 2 , in which the first term comes from power counting, 

that is the fluctuation level that would be expected from Kolmogorov scaling, while the second, 

_ i 

associated with production of intermittency, comes explicitly from the factor a m 2 in front of £, i.e. 
from the finiteness of the correlation time [12]. 

Coalescence processes 

In this case, the starting kinetic equation, can be written in the form, after discretization: 

^ x — l oo 

N(x,t + 1) = N(x,t) + - W(y, x-y-t)-^ W(x, y; t), (9) 

y=l y=l 

where W(x,y;t) = ■ w(x,y — ► x + y,t;j x ,j y ), with w equal to one, if particles j x and j y 
coalesce into an object of mass x + y, and to zero if they do not. The particles are supposed 
distributed at random over scales below that of the mean free path. The reason for this assumption 
is to focus on production of fluctuations at macroscopic scales, rather than from problems in the 
process of coarse graining. 

Dividing the volume V, into K slots of size much larger than that of the particles, the probabil- 
ity of having W2(x,y) slots with two particles, will be a multinomial distribution, that in the limit: 
K ^> N(x),N(y) ^> W2(x,y) >■ 1, will be approximated by a multivariate gaussian. Notice that 
the inequality K ^> N(x) is nothing else than the diluteness condition, while W2(x,y) ^ 1 is the 
large sample condition necessary for coarse graining. The probability of having W(x, y) coalescence 
events will be a gaussian as well in the hypotheses: 1 <C W(x,y) <C W2(x,y), so that, analogously 
to Eqn. (4), it is possible to write: 

V[W] * J] ex P ( - 2^ym^(^ w ^ w ^ (10) 

where: a(x,y) = a(x,y — ► x + y) = V~ 1 a(x,y)v(x,y), with a the coalescence cross section, and 
v an effective velocity; [Ka(x,y) is the coalescence probability for two particles in the same slot]. 
In the continuum limit, the average of Eqn. (9) gives immediately Eqn. (2), while the fluctuations 
are governed by the stochastic equation: 

/*oo 

d t n(x,t) + n(x,t) / dya(x, y)N(y) - £(x, t) 
Jo 



dya(x - y, y)N(x - y)n(y, t) - N(x) / dya(x, y)n(y, t), (11 



where: 



o Jo 



oo 



£(x,t) = l - I dyAW(x -y,y;t)- / dyAW(x, y; t). (12) 
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Focusing again on power law kernels: a(ax,ay) = a^a(x,y), with: a(x, y)\ x ^>y — x s y@~ s , there 
will be a constant current solution: N(x) = x~"< with 7 = ^-j^, which will be valid if the infinities 
in Eqn. (2) cancel one another, i.e. when: 8 < 

The noise correlation < £(x,t)£(x' ,t') > is obtained from Eqns. (10) and (12) and, as in the 
case of fragmentation, contains a singular piece proportional to 8(x — x'): 

< £(x,t)Z(x',t') > = 8{t - t')[E s (x)6(x - x') + E(x,x% (13) 

where: "E s (x) = y dya(x — y, y)N(x — y)N(y) + J °° dya(x, y)N(x)N(y) and "E(x, x') = a(x', x — 
x')N(x')N(x — x') for x > x' . This leads to the expression: 

3 

< n(x,t)n(x',t') > s = -8{x - x')N(x) exp(-T(a;)|i - t'\), (14) 

where: T(x) = J °° Aya(x,y)N(y). There are now two possibilities: 
A: if 8 < ^y^-, the correlation time T(x) is finite, and the same situation of fragmentation, with: 

arises. Systems in which particle interactions are dominated by fluid 
mechanics are likely to behave like this [13]. 
B: if > 8 > ^-5^, the correlation time decreases with x, like a power of xo/x, with xq 

the smallest particles in the systems. The final outcome, at large x/xo, is one of singular 
fluctuations behaving like a white noise in time, with vanishing amplitude ~ T(a;) _1 . The kernel 

— 1 / 3 -L_ — 1 / 3 

for Brownian coagulation: a(x,y) = ^1/3^1/3 [3,14], and the product kernel: a(x,y) = 

(xy) s , both lead to this behavior. 
Physically, the situation in case (B) is a consequence of the process being dominated by deposition 
of large numbers of the smallest particles. This leads to a steady growth in the mass of each 
agglomerate, with fluctuations vanishing, as the size of the deposited particles and their number, 
go respectively to zero and to infinity. 

The non-singular part of the fluctuations becomes important if one wants to verify the validity 
of Eqn. (2). In case (B), it can be shown that the non-singular correlations: < n(x,t)n(x' , /') > NS 
are continuous in x and x', so that, if 8 < "M^, the divergences at y = in the integrals in 
Eqn. (11) cancel one another. In case (A), problems with divergences are absent to start with. 
Hence, power counting can be used in both cases without any problem, and: < n(x,t)n(x' ,t) > N s~ 
x~ x N(x)g(x I x'). The condition for correlations being negligible: < ^\^^^ as i ^ 00, 
becomes therefore: N(x) ~ a; -7 , with 7 < 1, i.e.: AN(x — ► 00) = 00 [15]. This is nothing else than 
the requirement of preservation of diluteness and large sample size as x — ► 00. 

To conclude: fluctuations in both fragmentation and coalescence processes satisfy the rule 
< An(x) 2 >~ AN(x) characteristic of uncorrelated processes. The origin of the non self- averaging 
properties of shattering systems, already observed in [7], and of agglomerating systems, for (3 > — 1 
[see definition after Eqn. (12)], appears to reside only in the decrease in the direction of the 

processes, of AN(x) ~ xN(x). Also the vanishing of < ^^ f ^ > , for fixed Ax and x — ► 00, in the 
case of deposition dominated coalescence, is a manifestation of the lack of particle correlations, and 
fluctuations are killed in this case simply by the large number of deposited particles. 

The first reason behind these results, is the independence of the individual breaking or coa- 
lescence events, and the absence of memory in the solutions of the kinetic equations. These are 
very common properties in the case of fragmentation. In the case of agglomeration, however, as- 
sumptions on the absence of ordered motions at small scales become necessary, which are natural 
in cases like e.g. that of Brownian coagulation, but are more difficult to justify in general. 
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The second reason for lack of particle correlation, is the Markovian nature of the processes 
considered. In the case of fragmentation, in particular, violations of the result: < An 2 >~ AN are 
obtained, as soon as the correlation time for the breaking events becomes non-zero. This effect is 
what lies behind the generation of intermittency in turbulence models of the type recently studied 
in [11-12]. 
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